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Summary:  The  purpose  of  the  paper  is  to  present  a brief  review  of  basic  theoretical  approaches  to  two- 
dimensional  (2D)  nonlinear  supercavitating  flows  in  the  framework  of  theory  of  jets  in  an  ideal  fluid.  In  this 
connection  discussed  are  Kirchhoff  and  Zhukovsky  methods,  Chaplygin  method  of  “singular  points”,  method 
of  integral  equation,  etc.  A simple  model  problem  of  a supercavitating  (SC)  flat  plate  at  zero  cavitation 
number  a = 0 is  chosen  to  illustrate  the  core  of  the  methods  and  their  comparative  effectiveness.  Some 
mathematical  aspects  of  open  and  closed  cavity  closure  schemes  are  studied  as  well  with  use  of  Chaplygin 
method  applied  to  a SC  plate  with  a spoiler  at  nonzero  cavitation  number.  An  influence  is  demonstrated  of  free 
and  solid  boundaries  onto  the  cavity  volume  and  hydrodynamic  characteristics  of  the  plate.  Mathematics  4.0 
software  is  used  as  a main  tool  for  the  flow  pattern  visualization  of  the  problems  under  consideration.  An 
analytical  exact  solution  is  presented  to  the  2D  nonlinear  flow  problem  of  an  arbitrary  supercavitating  foil  and 
numerical  results  are  discussed. 

1.  Basic  assumptions  of  theory  of  jets  in  an  ideal  fluid 

Theory  of  jets  in  an  ideal  fluid  is  appeared  to  be  one  of  the  best-studied  fields  of  theoretical  hydrodynamics 
dealing  with  flows  confined  by  free  and  solid  boundaries,  the  pressure  constancy  condition  being  satisfied  on 
the  former.  H.  Helmholtz  [8]  and  G.  Kirchhoff  [12]  were  the  first  to  formulate  and  solve  some  relatively  simple 
problems  of  the  theory  of  jets.  Nowadays  one  can  find  brilliant  surveys  of  the  advances  and  development  in  the 
theory  in  books  by  Birkhoff  & Zarantanello  [1],  Gilbarg  [5]  and  Gurevich  [7]  (see  also  English  translation  of  the 
book  [6]).  The  works  by  Tulin  [23],  Terentev  [21]  and  Maklakov  [16]  should  also  be  distinguished. 

As  it  follows  from  the  name  of  the  theory  itself,  we  suppose  a fluid  to  be  an  ideal  one.  For  simplicity  we  neglect 
the  gravity  influence  and  specify  that  the  flow  be  2D.  steady  and  incompressible.  The  flow  has  a velocity 
potential  p if 

v = grad  p . 

where  v = vx  + ivy  is  the  total  velocity  vector  in  the  flow.  The  harmonic  function  p(x,  y)  is  the  real  part  of  an 
analytical  function  of  complex  potential  (or  characteristic  function)  w{z)  = p + iip,  where  z = x + i y and  ( x , y) 
denote  axes  of  rectangular  Cartesian  coordinate  system.  The  imaginary  part  of  w is  called  stream  function  %b, 
the  velocity  vector  being  given  by 

dip  d'tp  dp  d'd 

Vx  dx  dy  ’ °v  dy  dx 

Function  has  a constant  value  on  each  separate  streamline  including  free  surfaces  while  p increases  along  a 
streamline  in  downstream  direction.  A so  called  conjugate  velocity  can  be  introduced  as  follows 

d'UJ  in 

— = vx  + ivy  = ve  , (1) 

which  is  a mirror  image  of  velocity  vector  v with  respect  to  y- axis.  In  formula  (1)  v denotes  absolute  value  of 
the  vector  v and  9 is  an  angle  made  by  the  vector  to  x-axis. 

As  a result,  a problem  in  question  is  considered  to  be  solved  when  the  complex  potential  function  w(z ) = w(x+iy) 
is  found.  As  is  customary,  011  all  the  solid  boundaries  of  a given  topography  the  kinematic  (flow  tangency) 
condition  is  applied,  which  requires  the  fluid  flow  to  be  tangent  to  the  surface  of  the  boundary,  that  is  derivative 
of  p with  respect  to  normal  vector  n be  equal  to  zero.  O11  the  other  hand  the  free  streamline  condition  of 
constant  speed  is  satisfied  on  all  the  free  surfaces.  The  main  difficulty  of  such  free  surface  flow  problems  is 
connected  to  nonlinearity  of  the  boundary  conditions  which  are,  moreover,  to  be  satisfied  on  the  surfaces  of 
unknown  geometry. 


Paper  presented  at  the  RTOAVT  Lecture  Series  on  '"Supercavitating  Flows”,  held  at  the  von  Kdrmdn 
Institute  (VKI)  in  Brussels,  Belgium,  12-16  February  2001,  and  published  in  RTO  EN-010. 
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2.  Model  problem  — cavitating  plate  at  zero  cavitation  number 

Consider  a 2D  cavitating  flat  plate  in  a uniform  flow  with  velocity  absolute  value  at  infinity  and  angle  of 
attack  a,  see  figure  1.  The  origin  of  the  Cartesian  coordinate  system  is  taken  at  the  plate’s  trailing  edge,  z-axis 
being  directed  downstream  and  y upwards.  There  is  an  incident  stream  with  speed  Woo  coming  from  the  left. 
The  region  occupied  by  the  fluid  is  bounded  by  the  plate  [AB\  and  by  two  semi-infinite  free  surfaces  (AC)  and 
(BC). 


C 


C 


Figure  1:  Flow  pattern  for  the  cavitating  flat  plate  at  zero  cavitation  number. 


Assume  the  velocity  on  the  streamlines  Vq  to  be  equal  to  Vao  and  therefore  the  cavitation  number 


P00  P() 

PV2cJ  2 


= 0. 


where  and  po  denote  the  pressure  at  infinity  and  within  the  cavity  correspondingly.  It  readily  follows  from 
the  Bernoulli  equation  for  an  upstream  reference  point  and  free  surface-boundary  point 

1 2 1 2 
Po  + ^’o  = Poo  + ' 

We  render  all  the  parameters  of  the  problem  nondimensional  by  a suitable  choice  of  scale  so  that  Voo  = 1 and 
plate  length  equals  to  unit  1 = 1. 


2.1.  Kirchhoff  method 


To  solve  the  problem  under  consideration,  it  is  convenient  to  determine  the  function 

CM  = Woo-y11 


(2) 


rather  than  the  complex  potential  w(z ) = tp  + bp  itself.  Indeed,  if  function  C,(w)  is  found  then  z(w)  is  easily 
defined  as 


/ 


Wood  z 1 

— - — aw  = — 
aw  Woo 


(3) 


Note  that  it’s  often  no  sense  in  reversing  function  z(w)  to  find  w(z),  all  the  more  so  that  the  operation  is  rather 
cumbersome. 

So  the  erux  of  the  problem  is  to  determine  w and  C as  functions  of  the  same  single  variable  so  that  equation  (3) 
can  be  integrated.  Since  function  C (w)  defines  a conformal  mapping  of  w-plane  onto  the  C-plane,  one  has  to 
determine  these  regions  and  map  one  to  another. 

If  one  choses  the  stream  function  i/;  to  have  the  value  zero  on  the  dividing  streamline  (CO AC)  and  (COBC)  and 
since  all  the  streamlines  are  by  definition  lines  of  constant  V;  and  therefore  become  lines  parallel  to  the  y-axis, 
then  the  iw-plane  has  the  form  shown  in  figure  2.  The  plate-cavity  combination  appears  as  a semi-infinite  slit 
in  this  plane.  The  velocity  potential  (p  at  the  edges  of  the  plate  (points  A and  B)  has  the  values  ipA  and  ips- 
It  is  easy  to  see  that 

C=^eM 

V 
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Figure  2:  The  complex  velocity  potential  w. 


Oil  [OA]  one  has  6 = — a and  Voo/v  ratio  varies  from  infinity  at  point  O to  unit  value  at  point  A.  On  the 
(■-plane  interval  [OA]  corresponds  to  a semi-infinite  line  inclined  at  angle  —a  to  the  positive  direction  of  x-axis 
and  starting  from  the  point  (^  = e-,“.  A symmetrical  line  inclined  at  angle  7r  — a and  starting  from  the  point 
(b  = — e-la  corresponds  to  another  part  of  the  plate  [ BO} . Angle  9 varies  from  —a  to  zero  and  ratio  Voo/v  = 1 
on  streamline  (AC).  Further  on.  on  streamline  ( CA ) angle  9 increases  from  zero  to  tt  — a and  Voo/v  remains 
the  same.  Thus  free  surface  on  the  z-plane  becomes  a half-circle  of  unit  radius  on  the  (-plane.  The  (-plane  is 
depicted  in  figure  3.  A stagnation  point  C on  the  z-plane  corresponds  to  (c  = 1. 


Figure  3:  The  ( plane. 


To  obtain  a solution  to  the  problem  under  consideration  it  is  sufficient  to  transform  w-plane  onto  (-plane,  both 
being  given  in  figures  2 and  3.  Such  a problem  is  fairly  simple  and  much  easier  than  that  of  determining  harmonic 
function  ip  in  the  region  z with  unknown  free  boundaries.  The  substitution  of  the  complicated  boundary  value 
problem  by  a simple  problem  of  conformal  mapping  is  the  core  of  the  Kirchhoff  method. 

The  appropriate  w — ( mapping  can  be  found  step  by  step. 

First,  one  transforms  (-plane  to  the  first  quadrant  of  the  auxiliary  r-plane,  see  figure  4,  by  conformal  mapping 


r = a 


c-u 

( + Ca 


( — e-i“ 

° ( + e-i“  ' 


B 

i C 

AO  B 

0 a 


Figure  4:  The  r plane. 

With  the  correspondence  between  ( and  r planes  we  have  for  point  C 

1 — 0 — (X 

i = a — . therefore  a = cot  — . 

1 + e_1“  ' 2 

Second,  a new  variable  T\  = r2  varies  in  the  upper  semi-plane,  see  figure  5. 
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Figure  6:  The  t plane. 
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As  a result  one  obtains 


o • 4 « 

2 sin  — 


ipB  = l V0 


a a 
2 + 7T  sin  — cos  — 


(5) 


The  total  force  F acting  on  the  cavitating  flat  plate  is  calculated  by  integration  of  pressure  distribution  coefficient 

Cv 


F = — 


1 pvt, 


ZB 

/ 


Cpdz , 


(6) 


where 


is  given  by  Bernoulli  equation 


n _ P-Poo 
p pvlo/2 


P / 2 2\  f-'^oo  I i 

P - Poo  = X («oo  -V  )=  — l-TTTo 


PK 


2 ' 2 
Taking  into  account  that  (dw  = Wood 2 on  the  plate,  we  find 


ICI5 


F = 


ip't'oo 


1 


1 

w 


£ dw . 


Using  the  same  technique  as  above,  we  arrive  at  the  following  expression  for  the  force  coefficient 


Cf  = Cd  + i Cl  — 


F 

PF 11/2 


ei(T/2-«)  27rshlQ 

4 + 7r  sin  a 


(7) 


It  is  obvious  from  what  was  dealt  with  above  that  Kirchhoff  method  is  quite  complicated  even  for  such  a simple 
problem  of  a cavitating  flat  plate  at  zero  cavitation  number.  Nevertheless  it  was  the  first  to  enable  one  to  solve 
free  surface  flow  problems.  A new  important  step  in  this  direction  was  made  by  N.  Zhukovsky  in  1890  [26].  He 
proposed  a new  approach  to  solution  of  theory  of  jets  in  an  ideal  fluid  problem,  which  was  significantly  improved 
as  compared  to  Kirchhoff  one. 


2.2.  Zhukovsky  method 

A method  proposed  by  Zhukovsky  can  be  applied  to  the  2D  free  surface  flow  problems  with  1-connected  flow 
region  bounded  by  only  straight  solid  boundaries.  A new  function  (so-called  Zhukovsky  function) 


u = log  C = — log 


dw 

Wood  2 


log  1-  id 

Voo 


(8) 


is  introduced  instead  of  Kirchhoff’s  ( = Voodz/dw  function.  Another  improvement,  following  to  Zhukovsky 
method,  is  that  one  has  to  connect  functions  w and  u through  variable  t varying  in  the  upper  semi-plane  rather 
than  directly  determine  a relationship  between  the  two  functions.  Being  aware  of  dependences  w(t)  and  u>(t), 
one  can  eliminate  variable  t and  hence  get  a Kirchhoff’s  solution  to  the  problem.  However,  this  operation  is 
quite  aimless  because  of  the  following  relationships 


enabling  one  to  solve  the  problem. 

The  meaning  of  a novel  approach  is  as  follows.  It  was  shown  in  the  previous  section  that  w-plane  is  bounded 
by  horizontal  lines  ip  = const  which  are  images  of  streamlines  in  the  physical  2-plane.  On  the  other  hand, 
w-plane  is  bounded  by  vertical  and  horizontal  lines  in  the  case  of  straight  solid  walls  in  2-plane.  Indeed,  such 
solid  boundaries  are  lines  of  constant  9,  that  is  imaginary  part  of  w,  while  free  surfaces,  where  pressure  (speed) 
constancy  condition  is  satisfied,  are  lines  of  constant  real  part  of  w.  Therefore  w and  w planes  are  bounded  by 
polygonal  segments.  As  a result,  both  w and  w domains  can  be  transformed  onto  the  upper  half  of  auxiliary 
t-plane  by  using  Schwarz-Christoffel  transformation. 
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Schwarz-Christoffel  formula  allows  one  to  transform  half-plane  1 into  interiority  of  a polygon  with  n vertexes 
(n-gon)  on  the  w-plane.  It  is  assumed  that  each  vertex  is  of  angle  a,  < 2ir,  i = 1, . . . , n and 

n 

= 7T (n  - 2) . 

i=  1 


It  this  case  the  formula  is  as  follows: 


C, 


n 

I n (t-U)a,/*-1dt  + C2, 


(10) 


where  points  ti  correspond  to  the  polygon  vertexes.  If  1*  = oo  then  term  (1  — tfe)a'e/7r_1  in  the  integrand  should 
be  excluded.  Right-hand  side  of  formula  (10)  includes  2 n + 3 unknowns,  namely,  n real  parameters  ti, . ..  , 1„, 
n — 1 of  n real  parameters  a,,  i = 1 , , n and  four  real  parameters  in  two  complex  ones  C\  and  C 2 . An  n-gon  is 
totally  described  by  2 n coordinates  of  its  vertexes  on  w-plane  what  generates  2 n conditions.  Therefore  one  can 
arbitrarily  chose  three  of  2n  + 3 unknown  parameters,  say  ti,  t2  and  <3.  This  fact  is  a version  of  a well-known 
postulate  that  any  conformal  mapping  is  determined  by  three  preassigned  boundary  points. 

Consider  the  same  model  problem  of  eavitating  flat  plate  at  zero  cavitation  number  formulated  above. 

We  already  have  got  the  relationship  between  w and  1 planes,  see  figures  2 and  6: 

w - <pBt2  . (11) 


It  is  to  be  emphasized  that  in  a neighborhood  of  point  1 = 0 

dw  „ , 

— = 2 <pBt  = 0(1) , 


where  0(1)  denotes  a quantity  which  is  of  the  same  order  of  magnitude  as  1. 

On  streamline  (AC),  v = v0 0 wherefrom  log(wlDC/e)  = 0 while  angle  9 increases  from  —a  up  to  value  zero  and 
therefore  u = — log('y/,y00)  + iff  varies  from  — i a to  0.  On  streamline  ( BC ) function  u>  varies  from  i(7r  — o)  to 
zero.  On  the  interval  [OA]  of  the  eavitating  plate  imaginary  part  of  u>  is  a constant:  Imw  = —a  while  its  real 
part  decreases  from  infinity  down  to  zero.  At  point  C angle  9 varies  in  a jump-like  manner  from  — a up  to 
7r  — a,  see  figure  7. 


Figure  7:  The  w plane. 


Thus,  domain  on  w-plane  corresponding  to  a region  occupied  by  the  fluid  on  2-plane  is  a triangle  OAB  with 
zero  angle  at  infinite  point  O and  right  angles  at  points  A and  B.  Assume  that  to  = 0,  1b  = 1 and  tc  = oo, 
see  figure  6,  (a  postulate  about  three  preassigned  boundary  points).  Then  Schwarz-Christoffel  formula  (10) 
reduces  to 

t 

u>(t)  = Ci  f df  = + C2 . 

J l)(l  + 0,2) 

where  o2  = \J<pa/<Pb  ■ At,  point  B,  w(l)  = i(7r  — a),  that  is  why  C2  = i(7r  — a).  The  integral  above  yields 


,,,  „ , 1 . (1  — a2)t  + 2a2  7r  \ 

W(l)  = Cl  ( --  aresm  +^)+^~  «) 
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At  point  A,  u)(—a2 


-i a,  that  is  why  C\  = — ia  and,  finally, 


u\  ■ • (!  - a2)t  + 2a? 

ui{t)  = 1 arcsm — — — (-  x ( — 


(12) 


(1  + a2)  t 

Using  a relationship  between  the  arc  sine  function  and  natural  logarithm,  expression  (12)  can  be  re-written  as 


(1  + a2) t 


+ i(--« 


ui(t)  - log 

ia2(t  — 2)  — it  + — l)(t  + a2) 

The  system  of  two  nonlinear  conditions 

oj(oo)  = 0 ; and  zb  = 

allows  the  two  unknowns  and  to  be  derived.  The  former  gives 

a[pa  , a 
a = .7 - cot  — . 

V <PB  2 

while  the  latter  merit  some  additional  explanation.  Equations  (9),  (11)  and  (13)  can  be  combined  to  give 


(13) 


ZB 


- / 


Wood  z dw 
dw  dt 


dt . 


wherefrom 


; - 

Voo 


O+rH 


. — 4 


a 
2 ’ 


what  coincides  with  expression  (5). 

The  total  force  F acting  on  the  cavitating  flat  plate  is  calculated  in  a manner  of  the  previous  section 


F = 


i Pvt, 


zb 

I 


Cp  dz  = — 


i f»’a 


1 

/ 


a 


dz  dw 
dw  dt 


dt. 


where 

Up  = 1-1  e^2^(f)  and  ^ = — ew(t> 

I I dw  Voo 

and  give  the  same  result  just  as  above. 

Nowadays  some  modifications  of  Zhukovsky  method  (so-called  hodograph  method)  is  widely  used  and,  for 
instance,  was  applied  to  flow  domains  bounded  by  not  only  polygonal  segments  but  by  smooth  curves  as 
well  [11].  Note  that  the  first  adaptation  of  the  classical  Zhukovsky  (hodograph)  method  to  flow  past  a body 
with  a curved  topography  was  considered  by  Levi-Civita  [15]  and  Villat  [24]. 


2.3.  Mixed  boundary  value  problem  method 

This  method  seems  to  be  easy  to  understand  and  is  very  often  used.  The  solution  is  derived  in  two  steps. 
The  first  one  allows  us  to  find  complex  potential  w(C)  or  its  derivative  dw/d ( by  mapping  w-plane  onto  the 
preassigned  auxiliary  (-domain  (in  a manner  of  Zhukovsky  method).  In  the  case  of  a cavitating  plate,  using 
w-plane  (depicted  in  figure  2)  and  airxiliary  ( upper  half-plane,  see  figure  8,  we  get 


dw 

dC 


N C, 


(14) 


where  N denotes  a real  parameter  to  be  determined. 

The  second  step  is  formulation  of  a mixed  boundary  value  problem  for  Levi-Civita  [15]  function 


. , dw 

wLc  = i log r- 

VoodZ 


9 + i log = 9 + it  . 

Voo 


(15) 


It  is  obvious  that  Levi-Civita  function  is  a product  of  the  imaginary  unit  — i and  Zhukovsky  function  (8). 
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Figure  8:  The  ( plane. 


With  the  correspondence  between  physical  z and  auxiliary  ( planes  the  mixed  boundary  value  problem,  see 
figure  8,  is: 

r = 0 as  ( < — 1 and  £ > b 
( — a , as  — 1 < £ < 0 ; 

e=\ 

I 7 t — a,  as  0 < ( < 6 . 

where  6 is  a coordinate  (to  be  determined)  of  image  of  point  B in  (-plane. 

The  solution  to  the  mixed  boundary  value  problem  (Riemann-Hilbert  problem  [17])  is  obtained  via  formulae 
proposed  by  Keldysh  and  Sedov  [10].  The  formulae  are  extremely  useful  for  both  nonlinear  and  linearized 
cavitating  flow  problem  and  used  also  in  many  other  fields  of  fluid  mechanics  and  so  merit  some  additional 
explanation. 

Below  presented  are  Keldysh-Sedov  formulae  in  the  case  of  the  upper  semi-plane.  It  can  be  readily  re-written 
for  another  (-regions  like  circle,  etc. 

The  Riemann-Hilbert  problem  is  formulated  for  a function  f(z ) = u{z)  + \v{z)  which  is  holomorphic  in  upper 
semi-plane  z = x + i y,  its  real  part  u(z)  being  given  on  a set  of  segments  [a*  , bk],  k = 1 ,N  on  the  x-axis 
and  its  imaginary  part  v(z)  being  given  on  the  other  portion  of  a;-axis.  Following  to  [10],  three  types  of  solution 
to  the  problem  exist: 

* solution  unbounded  at  all  points  Ofc,  b * (oo  — oo  class): 


HO 


1 

7ri  R(z) 


Homo 

c-z 


d(  + i 


Homo 

C-z 


d(  + Pn{x) 


(16) 


where  a\  = ajv  + 1,  Pn('x)  is  a polynomial  of  degree  N or  N — 1 depending  on  behaviour  of  f(z)  at  infinity,  and 


mo  = 

* solution  bounded  at  all  points  at  and  unbounded  at  all  bt  (0  — oo  class)  as  /(oo)  = 0: 


\|  n(z~ak)(z-h)  ; 

\|  k=  1 


H0  = -^f\ 

7T1  Rb(z) 


' 1ST  hfc  /V  ak+i  ^ 

y[00mi0  iY  [viomio  l 

h C-zRaiO  h C-*R*(0  7 

k.  ak  Ok  / 


where 


Ra(0  = 


N 


w n (z~a0 > m(z) = 

\ k= i 


\ 


N 


Y[{z~l>k)  ; 


k=  1 


★ solution  bounded  at  all  points  at,  bk  (0  — 0 class)  as  /(oo)  = 0: 

bk 


f(0  = 


N k , . N ^ 

y f y r 

C-zR(0  H J 


HO  d( 


mo 

\^k[  t-zRw  ' mo 


(17) 


(18) 
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The  latter  solution  exists  if  and  only  if  the  relationship  holds 

bk 


w afc+l 

/'  v(0 


£/««+■£/ 


R(z) 


dc  = o . 


(19) 


Coming  back  to  solution  of  the  model  problem,  one  has  to  use  the  latter  formula  (18)  along  with  additional 
condition  (19)  to  derive  wlc(C)  for  the  free  surfaces  detaches  smoothly  from  the  edges  of  the  plate  A and  B and 
therefore  velocity  absolute  value  is  finite  there.  Thus,  using  boundary  conditions,  one  arrives  at  the  expression 


b 

"lc(C)^~V(C  + 1)(C-6)  / 

-1 


0(t)  cl  t 

y/(i  + t)(b-t) 


which  yields 


Condition  (19)  gives 


^lc(C)  = ilog 


Vb= ( - \/6y/r+C 

yjb  — C + y/b  v/T+X 


7 r .1—6 

— a H b arcsin - 

2 1 + 6 


— a . 


or  6 = tan2 (a/2). 

Expressions  (14)  and  (20)  give  analytical  solution  to  the  problem.  Conformal  mapping  is 


exp(iwLc(C))dC  • 


(20) 


(21) 


Condition  zb  = z(b)  = Ze1(-7r  enables  one  to  find  unknown  parameter  N in  the  form 

, t 8 / c x i (x 

N = cos  — . 

4 + 7r  sin  a 2 

Taking  into  account  that  ips  = w(b)  = Nb2 /2.  one  obtains  the  relationship  coinciding  with  (5). 

Note  that  the  method  of  mixed  boundary  value  problem  is  effectively  applicable  to  free  streamline  problems 
with  curved  topography  of  an  obstacle.  Actually  in  this  case  just  9(z ) is  given  and  function  9(()  = Rbwlc  is 
unknown  on  £-axis  because  it  is  not  a constant  and  the  relationship  z(()  is  yet  to  be  found.  A special  integral 
equation  is  formulated  is  a manner  presented  below  in  section  5. 


2.4.  Chaplygin  method  of  singular  points. 

Chaplygin  method  of  singular  points  based  on  the  idea  of  determining  of  a holomorphic  function  in  flow  domain 
and  therefore  in  auxiliary  complex  plane,  the  function’s  zeros  and  poles  being  known  and  Liuville’s  theorem 
applied.  The  theorem  guarantees  uniqueness  of  a solution  obtained.  The  auxiliary  region,  being  the  image  of 
the  flow  region,  is  chosen  to  be  bounded  by  straight  lines,  segments  and  circular  arcs  so  that  one  could  cover 
all  the  auxiliary  complex  plane  by  mirror  images  of  the  initial  domain  of  atixiliary  variable  range.  Quadrants, 
rectangles,  disk  sectors,  etc.  are  usually  used  as  such  ranges  depending  on  special  features  of  a problem  under 
consideration. 

To  get  the  desired  solution  to  the  problem,  one  has  to  determine  two  analytical  functions,  namely  derivatives 
of  the  complex  potential  with  respect  of  physical  z and  auxiliary  u variables  as  functions  of  u: 


dw 

(lw, 


(m)  and  xiu) 


dw 
Mood  z 


V 

Voo 


e 


-w 


Conformal  mapping  which  transform  flow  region  into  an  auxiliary  one,  then,  is 


fdz  f dz  dw 

ZW  = / = / 1 z~du. 

J au  J aw  au 
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Chaplygin  practically  introduced  some  basic  postulates  (though  he  had  never  published  them)  which  of  consid- 
erable assistance  while  finding  a derivative  of  w.  If  w(u)  transforms  w-region  into  w-region  then: 

• a third-order  pole  of  function  dw/du  corresponds  to  a jet  of  infinite  width; 

• a second-order  pole  of  function  dw/du  corresponds  to  a jet  of  semi-infinite  width  (that  is.  jet  bounded  by  a 
free  or  rigid  surface).  Such  a semi- infinite  jet  is  called  “ocean”  in  [1]; 

• a first  order  pole  of  function  dw/du  corresponds  to  a jet  of  a finite  width; 

• in  points  of  auxiliary  w-plane  corresponding  to  a vertex  of  solid  wall  bounding  flow  region,  function  x(y)  has 
an  exponential  zero  or  singularity  depending  on  angle  of  the  vertex; 

• function  dw/du  has  a zero,  more  precisely  behaves  like  0(t  — t0),  in  point  to  corresponding  to  a point  on 
physical  2-plane  where  a streamline  divides.  It  can  be  a stagnation  point  or  leading  edge  of  stagnation  zone, 
etc.; 

• conformality  condition  for  derivative  of  complex  potential  dw/du  often  does  not  satisfied  at  points  where  free 
streamlines  spring  from  the  body  (jet  detachment  points).  To  avoid  this,  one  has  to  chose  a form  of  auxiliary 
w-region  so  that  its  boundary  in  points  corresponding  to  jet  detachment  points  makes  right  angles. 

• function  dw/du  has  a zero  or  pole  (depending  on  specific  features  of  the  problem)  in  points  where  conformality 
condition  does  not  satisfied  for  w — u transformation. 

Now  apply  the  method  of  singular  point  to  the  model  problem  considered  above  with  Kirchhoff  and  Zhukovsky 
methods,  see  figure  1.  The  first  quadrant  of  auxiliary  w-plane  is  chosen  as  a region  to  be  transformed  into  the 
flow  domain,  see  figure  9.  As  is  customary,  the  coordinates  of  three  points  O,  B and  C on  the  boundary  of 
the  first  quadrant  of  w-plane  are  chosen  arbitrarily  and  location  of  point  A is  to  be  found.  Actually  the  figure 
coincides  with  figure  4 for  r-plane  used  in  section  ‘Kirchhoff  method’. 


i C 

_A O B 

0 a 

Figure  9:  The  auxiliary  w-plane. 


Following  to  the  main  postulate  of  Chaplygin  method,  find  both  derivatives  of  w with  respect  of  2 and  w as 
functions  of  u = { + i rj. 

Function  x(u)  = dw/(voodz)  is  a bounded  function  in  the  flow  region.  It  becomes  zero  in  just  one  stagnation 
point  O,  where  u = a.  On  the  plate  [AOB]  corresponding  to  real  {-axis  of  u 


argx(w)  = —0{u)  = 


a , as  0 < u < a ; 
a — 7r , as  a < u < oo , 


while  on  free  streamlines  (AC)  and  ( BC ) (imaginary  rj- axis  of  u),  |x(w)|  = 1.  Due  to  the  Schwarz  principle 
of  symmetry,  all  the  zeros  of  x(w)  becomes  zeros  in  corresponding  mirror  points  with  respect  to  the  imaginary 
axis  and  poles  in  corresponding  mirror  points  with  respect  to  the  real  axis.  That  is  why,  being  continued  onto 
the  whole  w-plane,  function  x(w)  has  a zero  at  point  w = a and  a pole  in  u = —a.  Then  the  following  ratio 


u — a 
u + a 


is  a bounded  and  holomorphic  function  everywhere  including  infinite  point  and  therefore  Is  a constant  due  to 
Liouville  theorem.  Thus 


x(«) 


d w 

Wood  2 


= Ni 


u — a 
u + a ’ 


where  N\  is  a constant  defined  by  condition  x(0)  = e‘“  which  yields  N\  = e'*“  7d  and  finally 


X(u) 


Wood2 


u + a 


(22) 
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Derivative  dw/du  has  a third-order  pole  at  infinite  point  C.  where  u = i,  and  zeros  at  stagnation  point  O 
(u  = a)  and  at  point  A (w  = 0)  where  conformality  condition  does  not  satisfied.  There  is  no  other  pole  or  zero 
in  the  first  quadrant  of  w-plane  for  dw/du.  On  continuation  of  the  function  into  the  whole  w-plane  by  mirror 
mapping  of  the  first  quadrant  with  respect  to  real  and  imaginary  axes  of  u (where  if’  = Im  w = const)  we  get 
additional  third-order  pole  at  point  u = — i and  zero  at  point  u = — a . Just  as  above,  the  ratio 

( u{u2  — a 2)  / dwA 

\ ('u2  + l)3  / dw  / 


is  a bounded  and  holomorphic  function  everywhere  including  infinite  point  and  therefore  Ls  a constant  due  to 
Liouville  theorem.  Therefore 


dw  u(u 2 — a2) 

dw  (u2  + l)3 


(23) 


where  N is  a real  constant. 

Two  expressions  (22)  and  (23)  allows  one  to  solve  the  problem  under  consideration  which  has  two  unknown  pa- 
rameters a and  N.  With  the  correspondence  between  z-jilane  and  w-plane.  the  conformal  mapping  transforming 
latter  into  former  is  as  follows 


s(«)  = fa Z-^L  Au  = f “(M  + a)2  dw  . 

{ 1 J dw  du  V ^ J (w2  + l)3 


(24) 


Note  that  (23)  and  (24)  can  be  readily  integrated  to  give 


w(u)  = 


N 


4(1  + a2)  \ w2  + 1 


o,2  _ „2 


(25) 


Two  conditions 


X(i)  = 1 • and  zB  = z(  oo)  = lei(”-a) 
generate  the  following  nonlinear  system  in  two  unknowns 


(26) 


i(a-ir) 


i — a N f u (w  + a) 


= 1; 


i + a ’ Woo  J (w2  + l)3 
o 


fu( 

J (u 


\rdu  = lei{7r-a) 


which  can  be  rewritten  to  yield 


a ,T  , . o a 16 

a = cot  - ; JV  = Woo  l sin  - — : 

2 2 4 + 7r  sm  a 


It  is  seen  from  (25)  that 


, , 1 .T  . o . . 1^T  o Ot  n (X 

ipB  = w(og)  = - N sm2  - , (pA  = w(0)  = - N cos2  - cot2  - , 


4 2 


4 2 2 


what  coincides  with  expression  (5). 

Taking  into  account  that  on  the  plate  [AB] 


(27) 


n i y2  i i i m2  i 1 dwdui 
c-  = 1 " - 1 “ Wo)l 


(28) 


the  total  force  F can  be  derived  (6): 


F = / Cpdz 


/- 


2 
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Figure  10:  Flow  pattern  for  the  cavitating  flat  plate  at  zero  cavitation  number  and  a = 50°. 


to  give  the  force  coefficient  in  the  form 


Cf  = 


CO 

if 


1 dw  dw  \ dz  dw 


du  = 


/ oo  oo 


dz  dz  ) dw  du 
o 

if  1 f dw  dw  _ \ 

~-i  Jisi;**)- 

Using  expressions  (22)  and  (23)  one  easily  obtains  coefficient  Cf,  see  (7). 


dw  dw 
dw.  dz 


du 


(29) 


Flow  pattern  z(u ) = x+iy  can  be  calculated  in  a parametric  form  from  expression  (24)  or  (26).  The  coordinates 
of  the  stagnation  point  O are  as  follows 


zo  = z(a)  = l e 


— 1 


sm  a ( a 

: 7T  — a + 2 cot  — + sm  2a 

4 + it  sm  a \ 2 


Flow  pattern  for  the  cavitating  plate  at  a = 50°  and  zero  cavitation  number  is  depicted  in  figure  10.  Stagnation 
point  O is  shown  as  well.  Figure  11  illustrates  behaviour  of  lift  and  drag  coefficients  Cl  and  Co  versus  angle 
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of  attack.  It  is  obvious  that  the  lift  equals  to  zero  for  a = 0 and  tt/2.  The  corresponding  curve  attains  its 
maximum  value  Cl  « 0.514  at  a*  « 39.26°.  Hydrodynamic  fineness 


* CD+Cf 

where  (7/  denotes  friction  coefficient  is  shown  in  figure  12  for  Cf  = 0.003  and  0.008.  The  former  curve  attains 
its  maximum  (11.237)  at  «maxi  = 2.5°  while  the  latter  one  (6.798)  at  amax 2 ~ 4.07°. 


Figure  12:  Hydrodynamic  fineness  versus  angle  of  attack  a for  friction  coefficient  Cf  = 0.003  and  0.008. 

It  is  to  be  underlined  that  Chaplygin  method  of  singular  point  appears  to  be  the  most  effective  while  solving  2D 
nonlinear  problems  of  theory  of  jets  in  an  ideal  fluid.  It  can  be  readily  applied  to  considerably  more  complicated 
problems  than  that  discussed  above  and  gives  quite  satisfactory  results. 


3.  A variety  of  cavity  closure  schemes 


Consider  a flat  plate  with  attached  cavity,  the  pressure  within  the  cavity  being  p = Po-  In  the  flow  region  the 
following  dynamic  condition  is  satisfied  p > po-  or  using  Bernoulli  integral 

v = Vo  on  the  cavity  boundary ; v < vq  in  the  ambient  fluid . 


Then  the  cavitation  number  a can  be  rewritten  as 


Poo  PO 


4--i>o, 


where  vo  is  an  absolute  velocity  value  on  the  cavity  boundary.  As  a result,  the  cavity  length  becomes  finite 
because  of  the  cavity  pressure  is  less  than  ambient,  including  that  at  infinity.  The  smaller  the  cavitation  number, 
the  larger  the  cavity  extent  and  in  a limiting  case,  as  a — » 0,  the  cavitation  flow  coincides  with  a streamline  one 
considered  above. 

The  Brillouin  paradox  [2]  is  well-known:  cavity  of  finite  length  with  closed  continuous  boundary  is  mathemati- 
cally impossible.  Indeed,  if  the  upper  and  lower  parts  of  the  cavity  make  a closed  contour  then  a stagnation  point 
should  appear.  That  is  impossible  due  to  the  requirement  of  a constant  absolute  velocity  value  on  the  cavity 
boundary.  That  is  why  a set  of  cavity  closure  models  has  been  developed,  for  the  definitions  and  properties  of 
which  reference  may  be  made  to  [1,  5,  7,  23].  Each  scheme  has  its  own  advantages  and  disadvantages  from  the 
mathematical  and  physical  viewpoint  and  should  be  chosen  depending  on  specific  features  of  the  problem  at 
hand. 

The  section  contains  a set  of  analytical  solutions  to  the  nonlinear  problem  of  a supercavitating  flow  past 
a flat  plate  of  chord  l,  with  angle  of  attack  a and  cavitation  number  a > 0,  the  cavity  closure  scheme  being 
varied.  The  schemes  considered  are  as  follows:  Tulin-Terentev  (single  spiral  vortex  termination),  Efros-Kreisel- 
Gilbarg  (re-entrant  jet  termination),  Riaboucliinsky  (symmetrical  plate  termination),  Tulin  (double  spiral  vortex 
termination),  Zliukovsky-Rosliko  (horizontal  plates  termination).  All  the  analytical  solution  were  derived  via 
Chaplygin  method  of  singular  points.  Numerical  results  were  obtained  in  Mathematica  for  Windows  computer 
mathematical  environment. 
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3.1.  Tulin— Terentev  scheme:  single  spiral  vortex  termination 

Consider  a cavitating  plate,  see  figure  13.  A cavity  closure  scheme  involving  single  spiral  vortex  is  accepted  as  a 
model  of  cavity  termination.  Tulin  [23]  was  the  first  to  propose  the  scheme  which  was  later  thoroughly  analysed 
by  Terentev  [7] . In  the  point  of  cavity  ‘collapse’  the  model  requires 


where  A > 0. 


Figure  13:  Flow  pattern  (0-plane)  and  auxiliary  w-plane  for  the  cavitating  flat  plate,  Tulin-Terentev  scheme. 


On  the  dividing  (‘zero’)  streamline,  where  ip  = 0,  the  following  relationship  holds  in  a neighborhood  of  point  C 


lim 


= ±00 , 


what  corresponds  to  two  spiral  streamlines  with  centers  at  points  C+  and  C . Sketch  of  the  flow  in  the  region 
of  the  cavity  termination  is  shown  in  figure  14. 


Figure  14:  Schematics  of  the  flow  in  the  vicinity 


Following  to  the  main  postulates  of  Chaplygin  method,  with  the  correspondence  between  physical  z-plane  and 
first  quadrant  of  auxiliary  m, -plane,  see  figure  13,  and  applying  Schwarz  principle  of  symmetry  to  cover  all  the 
w.-plane,  one  obtains  derivatives  of  complex  potential  w in  the  form 


dw 
d 0 


vo  e 


i(a  — 7r) 


u — a 
u-\-  a 


exp 


dw  u(u2  + l)(u2  — a2) 

d'u  (■ u 2 — u%P)2{it?  — Moo2)2  ’ 


(30) 
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wherefrom 


z{u) 


N_  i(7r_a)  f u(v-  + !)(><  +«)2 

v0  ./  («2-«L)2(«2-^2)2 

o 


x exp 


2 bu  \ 
u 2 + 1 ) 


d u . 


(31) 


It  is  easy  to  see  that  stagnation  point  O generates  zeros  at  point  u = a for  both  the  derivatives,  point  at 
infinity  D corresponding  to  u = mnfty  = c+ id  generates  second-order  pole  of  dw/du  and  at  point  u = i (image 
of  points  C,  C+  and  C-).  dw/du  = 0(u  — i)  and  uj  = 0(u  — i)-1. 

Three  conditions 


dw 

dz 


( Uqo  ) 


— ^OO  ? 


ZB  = z( oo)  = lei(7r-a) ; 


U oo 


0 


(32) 


give  a nonlinear  system  of  five  equations  in  five  unknown  parameters  of  the  problem  a,  b,  c,  d,  and  N to  be 
derived.  The  latter  condition  implies  an  existence  of  a contour  completely  surrounding  the  plate  and  the  cavity 
on  the  first  Riemann  sheet  and  actually  is  a cavity  closure  condition.  Terentev  was  the  first  to  propose  such 
a relationship  [20]  which  can  be  used  to  tell  so-called  closed  cavity  closure  schemes  from  ‘open’  ones.  All  the 
calculations  are  effectively  accomplished  with  use  of  Mathematics  for  Windows  software. 


Figure  15:  Flow  pattern  for  the  eavitating  flat  plate  at  a = 1 and  a = 30°. 

Force  coefficient  Cf  is  derived  in  a similar  manner  of  formula  (29): 

OO  

1 f dw  dw 

zb-  za o / I-  -r-du 

Vq  I du  dz 
o 

Figure  15  illustrates  the  flow  patter  for  a = 30°  and  a = 1.  It  is  to  be  emphasized  that  single  spirals  at  point  C 
shrink  very  rapidly  and  so  the  region  of  double-sheeted  flow  is  quite  small.  The  flow  region  in  the  vicinity  of 
point  C~  is  shown  in  figure  16.  Note  that  vertical  distance  S between  C+  and  C is  directly  connected  with 
drag  coefficient  Co- 

CD  = (l+a)j. 


CF  = -i 


. 1 + a 


3.2.  Efros— Kreisel— Gilbarg  model:  re-entrant  jet  termination 

The  flow  pattern  is  depicted  in  figure  17.  The  model  proposed  almost  simultaneously  by  Efros,  Kreisel  and 
Gilbarg  [1]  involves  a re-entrant  jet  of  unknown  width  8 and  direction  /i  at  infinity  point  C which  lies  on 
the  second  Riemann  sheet.  Gurevich  [7]  was  the  first  to  successfully  apply  the  scheme  to  the  problem  of  a 
supercavitating  plate  which  is  perpendicular  to  an  inflow.  An  additional  stagnation  point  E is  appeared  in  the 
vicinity  of  the  cavity  trailing  edge. 
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Figure  16:  Flow  region  in  the  vicinity  of  lower  single  spiral  vortex. 


Figure  17:  Flow  pattern  (0-plane)  and  auxiliary  u, -plane  for  the  cavitating  flat  plate,  Efros-Gilbarg  scheme. 


Adopting  first  quadrant  of  auxiliary  v-plane  as  an  image  of  the  flow  region,  with  the  correspondence  between 
physical  0 and  auxiliary  u planes  (see  figure  17)  solution  to  the  problem  is  as  follows: 

d w = ei(a_7r)  v - a (u  - m.0)(m  - W)  . 

d0  u + a (u  + uq)(u  + Mo)  ’ 

dw  _ u(u2  — a2)(u2  — Uq)(u2  — Mo2) 
dw  (m2  + 1)(m,2  - m20)2(m2  - m7^2)2 

where  u 0 = b + ic  and  Moo  = d + if-  These  two  expressions  can  be  combined  to  give 


0(M) 


U 

A [ Un-c)  f u{u  + a)2  (m,  + m0)2(m,  + Mp)2  , 

J e«2  + D (m2-«4)2(m2-m^2)2  " 


(35) 


The  same  conditions  (32)  yield  a nonlinear  system  of  five  equations  in  six  unknowns  a,  b,  c,  d,  f and  N.  It  is 
easy  to  see  that  the  solution  is  not  unique  for  the  number  of  unknowns  is  greater  than  the  number  of  conditions. 
Nevertheless,  one  can  use  an  additional  condition  connected  with  the  direction  of  the  re-entrant  jet  at  the  point 
of  infinity  C,  see  figure  18.  Both  the  asymptotic  analysis  and  numerical  results  for  the  cavitating  flat  plate  have 
shown  angle  fi  not  to  significantly  affect  the  hydrodynamic  coefficients  and  even  most  of  the  flow  pattern.  That 
is  why  it  seems  reasonable  to  consider  the  parameter  as  a given  one.  Thus,  an  additional  condition  is  as  follows 


dw 

d^ 


(i) 


Vo  e 1/1 . 


Using  residue  theory,  one  arrives  at  the  following  relationships  (which  are  correct  for  an  arbitrary  cavitating 
hydrofoil): 


2 g 


v0 


Cd  = 7 1 cos  u 

Mo J ' 


2g 


MO 


Cl  = — -7 sin  /i 

Moo  l 1 " 


(36) 
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where  q denotes  the  flow  rate  in  the  re-entrant  jet  C and  F is  circulation  along  a large  contour  completely 
surrounding  the  cavitating  foil  and  the  cavity  and  enclosing  most  of  the  flow.  Note  that 


res  = 


/ 


U-oo 


T + iq 


and.  moreover, 

(1  + a2)  ((1  + b2  — c2)2  + 4 b2c2) 

q = ttN ^ — 

2 ((1  + d2  - f2)2  + 4 d2f2)2 


Figure  18:  Flow  pattern  for  the  cavitating  flat  plate  at  a = 1,  a = 30°  and  /i  = 180°. 

Another  form  of  the  force  coefficient  is 

CF  = — (-yr+^'Ini  (res)  + ires)  . (37) 

^OO  * ' ' 

Flow  pattern  for  the  cavitating  flat  plate  at  a = 1 and  a = 30°  is  shown  in  figure  18  along  with  hydrodynamic 
coefficients  and  position  of  stagnation  points  O and  E.  The  direction  of  re-entrant  jet  is  chosen  to  be  /i  = 180°. 


3.3.  Riabouchinsky  model:  symmetrical  plate  termination 

The  model  introduced  by  D.  Riabouchinsky  [19]  assumes  the  flow  region  to  be  symmetric  and  so  the  cavity 
terminates  on  the  symmetrical  plate.  Another  version  of  the  scheme  was  proposed  [21]  with  central  symmetry, 
see  figure  19,  which  is  more  convenient  from  the  viewpoint  of  mathematical  analysis.  Adopting  this  scheme, 
transform  flow  region  into  the  rectangular  7r/2  x 7t|t|/2,  where  t is  purely  imaginary,  on  the  auxiliary  w-plane, 
see  figure  19. 


Figure  19:  Flow  pattern  (2-plane)  and  auxiliary  M-plane  for  the  cavitating  flat  plate,  Riabouchinsky  scheme. 
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Conjugate  velocity  dw/dz  has  constant  absolute  value  vq  of  both  vertical  sides  of  the  rectangular,  while  on 
the  horizontal  sides  its  argument  is  a step  function  with  a jump-like  behaviour  in  stagnation  points  O and  O' 
(argument  varies  by  7r).  That  is  why  dw/dz  has  zeros  at  points  u = a and  u = n — a + 7rr/2.  Using  Schwarz 
principle  of  symmetry,  one  arrives  at  the  conclusion  that  dw/dz  is  a doubly  periodic  elliptic  function  with 
periods  7r  and  7rr.  Basing  on  the  Liouville  theorem  and  accepting  notation  of  Whittaker  & Watson  [25]  one  can 
write  down  such  a functions  as  a ratio  of  elliptic  theta-function  i = 1, ....  4 which  nome  0 < q < 1 is  real. 
Function  dw/du  is  real  on  horizontal  sides  of  the  rectangular  and  must  be  purely  imaginary  on  vertical  sides. 
That  is  why  this  derivative  is  also  a doubly  periodic  elliptic  function  with  periods  n and  ttt.  The  point  D 
at  infinity  becomes  the  rectangular  center  u = 7r/4  + 7rr/4  where  dw/du  has  a second  order  pole.  Stagnation 
points  generates  simple  zeros  and  so  do  points  where  conformality  condition  does  not  satisfied,  that  is  all  the 
rectangle  vertexes. 


Figure  20:  Flow  pattern  for  the  cavitating  flat  plate  at  a = 1,  a = 30°. 

The  solution  to  the  problem  is  as  follows 

dw  Ua—ir)  ih{v  - a)ih(u  + a) 

— ft  ' ’ . 

dz  i)i(u  + a)$3(it  — a) 

(38) 

^ = N x Mu  ~ a)Mu  + a)Mu  - a)Mu  + a)  ■ 

wherefrom 

,(«)  = " =«*-)  / du . (39) 
Vo  J M 2m) 

o 

An  advantage  of  the  model  proposed  is  that  it  always  involves  less  number  of  unknown  parameters  then  the  other 
schemes.  Indeed,  one  has  to  determine  just  three  unknowns  a,  r and  N in  the  problem  under  consideration 
instead  of  5 — 6 in  other  cases.  Two  former  out  of  three  conditions  (32)  give  the  nonlinear  system  of  three 
equations  in  three  unknowns,  where  Moo  = 7t/4  + 7tt/4,  the  latter  condition  (32)  being  satisfied  due  to  the 
central  symmetry  of  the  flow  region. 

Force  coefficient  and  flow  pattern  are  calculated  using  expressions  (33)  and  (39)  correspondingly.  Figure  20 
illustrates  some  numerical  results  obtained  in  Mathematics  package  for  the  same  set  of  the  flow  parameters 
which  was  considered  in  the  previous  sections. 

3.4.  Tulin  model:  double  spiral  vortex  termination 

The  scheme  was  proposed  by  M.  Tulin  [23]  and  involves  two  double  spiral-like  streamlines  at  the  trailing  edge 
of  the  cavity.  Such  a double  spiral  vortex  flow  is  known  to  be  a unique  opportunity  to  conjugate  smoothly  two 
streamlines  with  different  absolute  velocity  values  [1],  see  figure  21.  The  scheme  is  appeared  to  be  especially 
efficient  for  confined  streamline  problems  with  one/two  free  surfaces  which  bounded  the  flow  region  [14,  22]. 
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Figure  21:  Flow  pattern  in  the  vicinity  of  a double  spiral  streamline,  from  [1], 


In  the  case  of  the  cavitating  flat  plate,  first  quadrant  of  auxiliary  w-plane  is  chosen  as  image  of  the  region 
occupied  by  the  fluid  in  2-plane,  see  figure  22.  The  wake  behind  the  cavity  begins  at  points  C and  E and  thins 
continuously  in  the  downstream  direction. 

Using  Chaplygin  method  of  singular  points  one  obtains  that  function  dw/dn  is  real  on  £-axis  and  purely 
imaginary  on  r/-axis.  It  has  zeros  at  stagnation  point  O (w  = a)  and  at  point  A where  conformality  condition 
does  not  satisfied  and  three-order  pole  at  point  at  infinity  D (u  = i). 


Figure  22:  Flow  pattern  (z-plane)  and  auxiliary  w-plane  for  the  cavitating  flat  plate,  Tulin  scheme. 


Another  derivative  dw/dz  has  a zero  at  stagnation  point  O and  behaves  in  a special  manner  at  points  C and  E. 
Consider  a function 


d wy 
d z J 


— e0+ilogu 


which  is  a step  function  on  £-axis.  Tracing  small  semi-circles  described  about  points  C (u  = ic  and  E (u  = id)  in 
the  counter-clockwise  direction,  one  arrives  at  the  conclusion  that  the  argument  of  the  function  has  an  increment 
(l/7r)(log(wo/w00)  and  (l/7r)(log(woo/wo)  correspondingly.  That  is  why 


where 


dz  u + a 


f u — ic  u + id 
\ u + ic  u — id, 


dw  u2  — a2) 

— = N — 

d u (u2  + l)3 


K = - log 

7T 


^oo 

^0 


(40) 


[ du  = 

A i(7T— a) 

f w(w  + a)2 

( u + ic  u — id  ^ 

J dw  dw 

1)0  J 

' (w2  + l)3 

\ u — ic  w + id,  J 

dw . 


and 
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The  problem  involves  four  unknown  parameters  a,  c,  d and  N.  On  the  other  hand,  just  two  ‘direct’  conditions 
can  be  imposed: 

“T~(^oo)  = Vqo  and  zb  — 2(00)  = le1^7V~cx^ . (41) 

d z 

Cavity  closure  condition  applied  in  the  previous  sections 


^du  = 0 
du 


does  not  satisfied. 

Thus  this  scheme  is  appeared  to  be  the  most  indeterminate  cavity  closure  model.  Nevertheless,  two  additional 
conditions  should  be  formulated: 

= 0 (42) 

?7=1 

to  close  the  problem.  The  former  implies  the  velocity  potential  to  be  the  same  at  the  end  of  the  cavity  at 
upper  and  lower  boundaries.  The  latter  specifies  the  angle  9 (direction  of  the  velocity  vector)  along  the  wake 
boundaries  be  minimal  at  point  D.  Note  that  both  conditions  are  arbitrarily  chosen  and  can  be  substituted  by 
a pair  of  others,  for  instance  yc  = Ve  and  d2w/d if  = 0 as  r/  = 1. 

Numerical  results  obtained  with  use  of  conditions  (41)  and  (42)  are  shown  in  figure  23. 
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Figure  23:  Flow  pattern  for  the  cavitating  flat  plate  at  a = 1,  a = 30°. 


3.5.  Zhukovsky— Eppler— Roshko  model:  horizontal  plates  termination 


The  model  is  often  called  ‘open’  one  and  assumes  that  free  streamlines  conjugate  smoothly  with  two  solid  semi- 
infinite plates  which  are  parallel  to  the  inflow.  The  absolute  velocity  value  monotonically  decreases  from  vq 
down  to  Voc  along  this  plate.  This  cavity  closure  model  is  very  effective  for  streamline  problems  either  under 
gravity  or  with  one/two  solid  boundaries  [7].  Zhukovsky  [26]  was  the  first  to  introduce  the  model. 

Following  to  Chaplygin  method,  transform  flow  region  into  the  rectangular  7t/2  x 7rjrj/2,  where  r is  purely 
imaginary,  on  the  auxiliary  u-plane,  see  figure  24.  Determining  zeros  and  poles,  one  arrives  at  the  following 
solution  to  the  problem 


dw  _ i(a_T)  gi (m-o)  . 
dz  0 f?!  (•«  + «)’ 

dw  r9i  (2m)  (?i  (m  — a)  $ \ ( u + a) 
du  $4  ( u — b)  $4  (m  + 6) 


(43) 


wherefrom 


z{u) 


N i(,r-a)  1 ih{2a)  d\{u  + a) 
M0  ,/  $!(m  - b)  #](m  + b) 


(44) 


10-21 


Figure  24:  Flow  pattern  (z-plane)  and  auxiliary  w-plane  for  the  cavitating  flat  plate,  Zlnikovsky-Roshko-Eppler 
scheme. 
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Figure  25:  Flow  pattern  for  the  cavitating  flat  plate  at  a = 1,  a = 30°. 


First  two  out  of  three  conditions  (32)  where  u0 0 = 6+  ttt/2,  along  with  two  additional  conditions 


din 

dz 


— Voo 


and 


(lloo)  = 0 


allow  four  unknowns  a,  b,  r and  N to  be  derived.  The  former  of  the  additional  conditions  implies  that  the 
semi-infinite  plates  are  horizontal  and  the  latter  specify  the  velocity  absolute  value  be  a monotone  decreasing 
function  having  its  minimum  at  point  D (u  = Moo). 

The  corresponding  flow  pattern  and  hydrodynamic  coefficients  are  shown  in  figure  25.  Note  that  this  open  cavity 
closure  scheme  give  significantly  smaller  cavity  length  than  so  called  closed  models,  such  as  Riabouchinsky, 
Tulin- Terentev,  Efros-Kreisel-Gilbarg.  As  a result,  the  hydrodynamic  coefficients  become  somewhat  larger. 


4.  Free  and  solid  boundaries 

All  the  previous  section  has  demonstrated  an  effectiveness  of  Chaplygin  method  of  singular  points  applied  to 
nonlinear  cavitating  problems  of  theory  of  jets  in  an  ideal  fluid.  This  section  contains  another  documentation 
of  the  fact  dealing  with  two  somewhat  more  complicated  problems:  first  one  being  that  of  the  flow  past  a 
cavitating  flat  plate  in  the  channel  and  second  one  - cavitating  plate  in  the  uniform  jet  of  finite  width.  It  was 
shown  by  Terentev  [21]  that  Zhukovsky-Roshko  and  Tulin  (double  spiral  vortex)  schemes  are  to  be  used  for  the 
problems.  We  shall  not  dwell  on  the  solution  procedure  and  just  describe  the  final  results  and  some  limiting 


cases. 
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4.1.  Cavitating  plate  in  a channel 

The  flow  pattern  z = x + i y and  the  auxiliary  zt-plane  are  depicted  in  figure  26.  Given  are  the  following 
parameters:  angle  of  attack  a,  plate  length  l,  channel  width  H,  cavitation  number  a and  a distance  between 
the  lower  wall  of  the  channel  and  point  A (trailing  edge  of  the  plate).  It  is  also  assumed  that  velocity  v\  at 
points  D i and  D:i  is  the  same  (though  it  is  an  unknown  value). 


Figure  26:  Flow  pattern  (z-plane)  and  auxiliary  u-plane  for  the  cavitating  flat  plate  in  the  channel. 
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where  Uoo  = 112  = d,2  + ttt/2,  m = di  + ttt/2,  i = 1,3. 

After  some  algebra  the  number  of  unknowns  (and  equations)  can  be  reduced  from  nine  down  to  six  and  even 
to  four,  but  the  system  remains  complicated. 

The  general  problem  has  lots  of  interesting  special  cases  as  hi  or  /13  tends  to  infinity  (one  solid  wall),  d\  = 0 
and  d‘i  = tt/2  (Kirchhoff  scheme  in  the  channel,  when  point  C coincides  with  D3  and  E with  D 1),  etc.,  see 
figure  27. 


4.1.  Cavitating  plate  in  a jet  of  finite  width 

The  flow  region  and  the  corresponding  auxiliary  w-plane  are  shown  in  figure  28.  Tulin  cavity  closure  model  with 
double  spiral  vortex  is  adopted  for  it  seems  to  be  the  most  efficient  from  the  mathematical  viewpoint.  The  jet 
width  is  assumed  to  be  If  and  submergence  of  the  plate  is  h.  Actually  submergence  is  a distance  between  a 
dividing  streamline  and  the  upper  jet  boundary  far  upstream.  The  velocity  on  the  jet  boundary  is  Woo. 


Figure  28:  Flow  pattern  (z-plane)  and  auxiliary  u -plane  for  the  cavitating  fiat  plate  in  the  jet  of  finite  width. 


If  one  choses  a first  quadrant  of  the  auxiliary  u- -plane  as  an  image  of  the  flow  region,  then  solution  to  the  problem 
is  found  through  Chaplygin  method  in  the  form: 


where,  as  beforehand 


and 


d w un-w\  u — a ( u — ic u + i d\K 

= VQe(  > — — — — ; 

d z u + a V u + ic  u — id.  J 


u(u2  - a2) 


— = N — 

dw  (u2  + 1 ) (u2  + e2)(u2  + /2) 


1 , Woo 

K — — log 

7T  V0 


. fdz  f dz  d w 

ziu)  = J S'1"',/  d^'d^d“- 


(48) 


Conditions 


dw 

dz 


(m  00)  = Woo!  Zb  = z(oo)  = le1^  ; arg  -j— (i)  = arg  (ie) : 

dz  dz 

w(i d)  = w(ic) ; ^ ^ (jjy  dw  = iHv^  ; ^ ‘f  = lhv°° 


allow  the  unknowns  a,  c,  d,  e,  / and  N/Voo  to  be  derived.  Again,  the  third  and  fourth  conditions  (which  specify 
the  velocity  vector  be  of  the  same  direction  at  points  E and  G at  infinity  downstream  and  complex  potential  be 
of  the  same  value  at  the  ends  of  the  cavity  C and  D correspondingly)  are  quite  artificial  and  can  be  substituted 
by  other  reasonable  relationships. 
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Figure  29:  Flow  pattern  for  the  cavitating  flat  plate  in  the  jet. 


Note  that  solution  (48)  produces  a solution  by  Larock  and  Street  [14]  for  the  cavitating  plate  beneath  the  free 
surface  as  its  special  case  as  H — ► oo  and  / — ► 1+.  Solution  for  the  cavitating  plate  in  unbounded  inflow  (40) 
is  also  a special  case  if  H , h — » oo  and  d , e — + 1+  as  well. 

Figure  29  illustrates  numerical  results  (flow  pattern  in  the  vicinity  of  the  plate  and  hydrodynamic  coefficients) 
obtained  with  use  of  Mathematics  for  the  cavitating  plate  with  a following  set  of  given  parameters:  H = 2, 
h = 0.75.  a = 30°,  a = 0.3. 

5.  Separated  free  streamline  flow  around  a body  with  a curved  boundary 

5.1.  Levi-Civita  approach 

Levi-Civita  [15]  was  the  first  to  propose  an  approach  to  solution  of  the  plane  free  surface  flow  problems  for 
obstacles  of  a curved  topography.  Such  problems  are  significantly  more  complicated  than  those  with  flow 
regions  bounded  by  polygonal  segments.  For  the  latter  the  imaginary  part  9(u)  of  the  Zhukovsky  function  oj(u) 
is  given  on  the  portion  of  the  auxiliary  w-domain  boundary  corresponding  to  the  straight  solid  walls  and  its  real 
part  (velocity  absolute  value)  is  given  on  the  other  portion  of  the  boundary.  In  the  case  of  a curved  obstacle 
all  information  one  possesses  is  just  the  function  6(z)  on  the  wetted  portion  on  the  body  (which  is  sometimes 
also  unknown)  and  velocity  absolute  value  on  free  surfaces.  The  main  difficulty  of  the  problem  is  that  it  is 
impossible  to  set  9(u)  without  conformal  mapping  z(u),  if  9(z)  is  not  a step  function. 

Following  to  Levi-Civita  approach,  flow  region  2 is  transformed  into  the  semi-circle  in  auxiliary  w, -plane  (one 
can  also  chose  another  domain  as  well)  and  Levi-Civita  function  (15)  is  written  down  in  the  form 

wlc  = wo  + - • 

where  u>q  denotes  Levi-Civita  function  for  a polygonal  body.  Additional  term  is 

OO 

n(«)=  Akuk> 

k——n 

where  A k are  real  coefficients  to  be  found  from  the  imposed  conditions  including  Brillouin  one  (smooth  detach- 
ment condition) . 

Brodetsky  [3]  calculated  a flow  around  a circular  and  elliptic  cylinder  for  a = 0,  three  coefficients  A ^ being 
taken  into  account.  I11  the  case  of  circular  cylinder  the  relationship  holds 

d 9 

ds  ~K’ 

where  s denotes  an  arc  coordinate  of  the  body  boundary  and  k is  its  curvature.  Collocation  method  can  be 
used  to  derive  the  coefficients  A fc.  Significantly  more  detailed  numerical  analysis  of  the  problem  was  done  in  [7] 
by  Terentev  and  Dmitrieva  for  an  arbitrary  cavitation  number  <7,  more  than  30  terms  being  taken  into  account. 
Two  points  were  determined  satisfying  the  Brillouin  condition.  The  points  correspond  to  two  maxima  on  curves 
Cx  versus  detachment  angle  7 for  a given  a,  see  figure  30. 
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Figure  30:  Drag  coefficient  of  a circular  cylinder  versus  detachment  angle  7,  from  [7]. 


5.2.  Method  of  integral  equation 

This  method  is  an  adaptation  of  mixed  boundary  value  problem  method  (see  section  2.3)  to  flow  around  an 
obstacle  of  a curved  topography.  By  a conformal  mapping  of  the  flow  region  into  the  semi-plane  (or  semi-circle, 
etc.),  the  solution  to  the  problem  is  reduced  to  an  integral  (or  integro-differential)  equation  of  Villat  type  [24] 
for  a function  X(u)  characterizing  the  angle  made  by  the  fluid  velocity  vector.  Using  another  form  of  such  an 
integral  equation,  Nekrasov  [18]  was  the  first  to  prove  both  existence  and  uniqueness  of  flows  around  a small 
circular  arc. 

Consider  the  method  applied  to  a free  streamline  problem  for  an  arched  contour,  see  figure  31,  and  transform 
the  region  in  2 plane  occupied  by  the  fluid  into  the  upper  semi-plane  £,  see  figure  8. 


^OO 


Figure  31:  Flow  pattern  for  the  cavitating  arc  contour  at  zero  cavitation  number. 


Then,  again. 


(49) 


where  N is  a real  parameter  to  be  determined  and  the  mixed  boundary  value  problem  for  Levi-Civita  func- 
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tion  (15)  arising  in  ((-plane: 


m = 


t = 0 as  £ < — 1 and  £ > 6 
— a(£) , as  — 1<£<0; 
7 t — a(£),  as  0 <£<6. 


where  6 is  a coordinate  (to  be  determined)  of  image  of  point  B in  ("-plane  and  «(()  is  the  body  angle  at 
point  ( = (.  relationship  a(s)  being  given,  where  s is  an  arc  coordinate  of  the  contour,  s = 0 at  point  A and 
s = S at  point  B. 

The  Keldysh-Sedov  formula  gives  the  solution  to  this  problem  in  the  0 — 0 class 


wlc(C)  = \/(C  + l)(C-&)  X 

7T 

which  exists  if  and  only  if 


i-i 


a{t) 


At 


+ 


o 

/ 


At 


\/ (1  +t)(b-t)  t-  ( J sj{l  +t){b-t)  t - C 


-a 


a(t)dt 


V(1  +t)(b-t)  2 


7T  .1-6 

+ — — arcsm  = 0 . 


1+6 


It  is  seen  that 


="<:*“(£ae) 


This  equation  can  be  used  in  conjunction  with  the  equality 


^ = exp(ia) 


to  give 


2 C 

f N ( 

s(£)  = / exp(— ia)dz  = — / exp { — ia(t)}  exp {iwLc(i)} t dt . 
J I'OO  J 


-1 

Substituting  the  solution  for  wlc  into  this  equation  finally  gives 

« 


N r 

s(0  = — G{t, b)  At, 

Vco  J 


(*) 


(50) 


where 


am  = + w -J)  +«.<■  - 1)  x „„  ( v-c  + »)(+  ‘>  , "W  <*  ] 

i+4  | * J v'a  + <)d-*) 1 - ( J 

This  integral  equation  along  with  relation  (*)  and  an  obvious  condition 

s(6)  = S 

allows  unknown  parameters  6.  N and  function  «(()  to  be  determined. 

5.3.  Arbitrary  2D  supercavitating  hydrofoil:  analytical  solution 

Both  the  Levi-Civita  and  integral  equation  methods  are  known  to  connect  with  significant  difficulties  while 
considering  a flow  past  an  obstacle  with  a boundary  of  large  curvature  k.  At  the  same  time,  hydrofoils  have 
such  a region  in  the  vicinity  of  the  leading  edge.  One  can  find  some  theoretical  analysis  of  the  nonlinear  problems 
of  an  arbitrary  supercavitating  hydrofoil  and  numerical  results  as  well  in  [9,  4]. 

That  is  why  an  approach  proposed  by  Maklakov  [16]  for  arbitrary  cavitating  hydrofoil,  see  figure  32,  Is  of  interest 
as  an  improvement  of  Levi-Civita  method. 
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Figure  32:  Flow  region  around  a cavitating  hydrofoil  (Tulin-Terentev  cavity  closure  scheme). 


Figure  33:  Auxiliary  (-plane  for  cavitating  hydrofoil. 


Consider  a cavitating  hydrofoil  which  boundary  is  given  by  relation  a = F(s),  where  s - arc  coordinate  (s  = 0 
at  point  A)  and  a is  the  tangential  angle  to  the  foil.  Tulin-Terentev  cavity  closure  scheme  with  a single  spiral 
vortex  model  of  cavity  termination  is  adopted  (subsection  3.1).  Following  to  Maklakov,  transform  flow  region  2 
into  the  upper  semi-circle  |t|  < 1,  Im  > 0 on  auxiliary  plane  t = £ + iry,  see  figure  33.  Stagnation  point  O 
corresponds  to  to  = exp(i<50),  point  C - to  the  origin  of  the  coordinate  system  in  t- plane  and  point  at  infinity  D 
has  the  image  at  too  = exp(iJoo)- 
Chaplygin  method  of  singular  points  gives  the  derivative 


dw 

d t 


kv0  f(t) 


kv  0 


t(t2  — l)(t2  — 2f  cos  A)  + 1) 


(t-tc 


>)  (t  - 1/too)  (t  - 1/too) 


(51) 


where  k is  a real  parameter  to  be  determined. 

Levi-Civita  function  is 

t2  — 1 

WLc(t)  = wo(t)  + M — f O(t) , (52) 

where  M denotes  a real  constant,  u>o(t)  is  Levi-Civita  function  for  a flat  plate  at  zero  cavitation  number 


W0(t)  = «0  - 7T  + i log 


t — exp(i<fo) 

1 — t exp(i<5o) 


and  Q(t)  is  analytic  and  continuous  function  in  upper  semi-circle.  That  is  why  conformal  mapping  can  be 
written  as 

t t 

z(t)  = / “j“  “jj  d*  = f exP  ji  vo(t)  + i M 1 ^ 1 + iO(t)  | f(t)  d t . 

1 1 

Levi-Civita  method  requires  the  expansion 


fi(t)  = ]TAfctfc, 

k=l 
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where  A k are  real  coefficients  to  be  found.  Maklakov  [16]  proposed  another  representation  of  0.  If 

tt(eis)  = \(6)+in(5) 

then  Schwarz  formula  for  a circle  gives  a relation  between  Q and  its  real  part  A given  on  the  semi-circle  boundary 


m-^j 


\(S)dS 


1 — 2 1 cos  8 + t2 


Imaginary  part  of  Q(t)  is 


7 r 

l4~f)  = c A = j {\{e)  - A(7))  fcot 


e+7  4£-7 

— cot  — - — 


de . 


Figure  34:  Zhukovsky  hydrofoil  (a)  and  relation  a = F(s). 

Using  expression  ds/d<5  = |dz/dt|  and  conformal  mapping  z(t)  write  down  the  relation  between  arc  coordinate  s 
and  polar  angle  8: 


s(8)  = k j g( 7)  exp(— CX  — 2M  SU17)  d7  , 


where 


g{ 7)  = 4 sin 7 (1  - cos (7  + <50))  x 


(1  + - 2 Too  cos (7  - (Soo))2  (1  + - 2-roc,  cos(7  + 4o))2 


It  is  easy  to  see  that  A = a — «o  and  therefore,  finally,  one  arrives  at  the  following  integral  equation 


A (5)  = — «o  + F J 3(7)  exp(— CX  — 2M  sin  7^7^ 


(53) 


Five  unknown  parameters  k.  M,  A0,  roo-,  800  are  derived  from  the  conditions 

77  dt  = 0 , ^-{too)  = ^00  , s(-l)  = S, 

at  d z 
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1 1.2  1.4  S 

Figure  35:  Drag  coefficient  Co  versus  wetted  length  S. 


1 1.1  1.2  1.3  1.4  S 


Figure  36:  Lift  coefficient  Cl  versus  wetted  length  S. 


where  S denotes  the  whole  (given)  wetted  length  of  the  hydrofoil.  The  latter  condition  can  be  substituted  by 
Brillouin  condition  at  point  B which  specify  the  curvature  of  the  cavity  and  hydrofoil  be  the  same  at  detachment 
point  B (t  = —1).  It  was  demonstrated  [1]  that  this  condition  is  equivalent  to  fi'(— 1)  = 0 which  yields  [16] 


S0 

7 r tan  — + 


71 

/ A'(7) 


tan  — d7  = 2-kM  . 


Naturally,  in  this  case  the  wetted  length  S is  unknown  because  ^-coordinates  of  point  B are  unknown  as  well. 
Some  numerical  results  for  a supercavitating  Zhukovsky  hydrofoil  are  shown  in  figures  34-36.  Zhukovsky 
hydrofoil  of  unit  chord  length  is  chosen  to  have  maximum  thickness  11.82%,  perimeter  P = 2.053  and  curvature 
radius  at  the  leading  edge  p = 0.016,  see  figure  34,  where  the  foil  geometry  and  relation  a = F(s)  for  the 
tangential  angle  to  the  foil  versus  arc  coordinate  s are  shown.  Figures  35  and  36  demonstrate  drag  and  lift 
coefficient  versus  wetted  length  S (Brilllouin  condition  is  not  imposed)  in  the  case  of  angle  of  attack  a = 5°  for 
a set  of  cavitation  number  a = 0,  0.1,  0.15,  0.3  and  0.5  (curves  T’  corresponds  to  a = 0).  It  was  found  that 
drag  (lift)  coefficient  attains  its  maximum  (minimum)  value  at  points  where  Brilllouin  condition  was  satisfied. 
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Conclusion 

A brief  review  of  the  main  theoretical  approaches  to  problems  of  the  theory  of  jets  in  an  ideal  fluid  was  conveyed, 
including  hodograph  (Kirchhoff  and  Zhukovsky)  methods,  Chaplygin  method  of  singular  points,  mixed  boundary 
value  problem  method,  some  modifications  of  Levi-Civita  approach  and  method  of  integral  equations.  Every 
analytical  solution  is  illustrated  by  numerical  results,  including  flow  patterns. 
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